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^ i Experiments in coherent spectroscopy correspond to control of quantum mechanical ensembles guid- 

ing them from initial to final target states by unitary transformations. The control inputs (pulse 
sequences) that accomplish these unitary transformations should take as little time as possible so 
as to minimize the effects of relaxation and to optimize the sensitivity of the experiments. Here, we 
present a novel approach for efficient control of dynamics in spin chains of arbitrary length. The 
approach relies on creating certain three spin encoded states, which can be efficiently transferred 
through a spin chain. The methods presented are expected to find applications in control of spin 
dynamics in coherent spectroscopy and quantum information processing. 

c ■ 

qh 1 Introduction 

> : 

According to the postulates of quantum mechanics, the evolution of the state of a closed quantum 
system is unitary and is governed by the time-dependent Schrodinger equation. This evolution can 
be controlled by systematically changing the Hamiltonian of the system. The control of quantum 
systems has important applications in physics and chemistry [1-4]. In particular, the ability to steer 
the state of a quantum system (or of an ensemble of quantum systems) from a given initial state to a 
desired target state forms the basis of spectroscopic techniques such as nuclear magnetic resonance 
(NMR) and electron spin resonance (ESR) spectroscopy ||] , laser coherent control [f| and quantum 
computing j^, ||. Achieving a desired unitary evolution in a quantum system in minimum time is 
an important practical problem in coherent spectroscopy. Developing short pulse sequences (control 
laws) which produce a desired unitary evolution has been a major thrust in NMR spectroscopy Qj. 
For example, in the NMR spectroscopy of proteins, the transfer of coherence along spin chains is an 
essential step in a large number of key experiments. 
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A spectroscopist has at his disposal a limited set of control Hamiltonians {Hj} (produced by 
external electromagnetic fields) that can be turned on and off to modify the net Hamiltonian of 
the system. There is a natural coupling (interaction) between the spins and in the absence of any 
external control Hamiltonians, the state of spin system evolves under this interaction or coupling 
Hamiltonian 7i c - The task of the pulse designer is to find the right sequence of external pulses 
interspersed with evolution of the system under the coupling Hamiltonian Ti c for different time 
periods, in order to create a net evolution or unitary transformation that transforms the state of the 
system from some initial to a desired final state in minimum possible time. 

Even for two coupled spins 1/2, the time-optimal transfer of polarization or of coherence is not trivial 
J?], ||] . Numerous approaches have been proposed and are currently used to transfer polarization 
or coherence through chains of coupled spins. Examples are the design of radio-frequency (RE) pulse 
trains that create an effective Hamiltonian [[l], [Hj, which makes it possible to propagate spin waves 
in such chains jl^, [l3| |TJ, |l5|, [l6|. In order to achieve the maximum possible transfer amplitude, 
many other approaches, that rely either on a series of selective transfer steps between adjacent spins 
or on concatenations of two such selective transfer steps [ll[] have been developed. 

Here, we consider a novel approach to control the transfer of coherence in spin chains of arbitrary 
length. The approach relies on the creation of a three spin encoded state and efficient propagation 
of this encoded state through the spin chain. Our method is based on variational ideas as captured 
by the theory of optimal control ||- In the present context of control of nuclear spin ensembles, 
we are interested in finding a sequence of RF pulses that will efficiently transfer any state (known 
or unknown) of a given spin in a spin chain to a desired target spin. Compared to conventional 
experiments, this new approach makes it possible to speed up the transfer rate by up to a factor of 
three, which suggests applications in NMR spectroscopy and experimental quantum computation. 

Without loss of generality, here we consider the problem of transferring the coherence of a spin at 
one end of a spin chain (label spin 1) to a spin at the opposite end of the chain (label spin n) in a 
spin ensemble. Consider an initial density operator po representing coherence or polarization on the 
first spin, which has the general form po = f 1 + ailix + ^Iiy + azl\z- The goal is to transfer this 
density operator to the operator ~1 + ciil nx + a 2 I ny + a^Inz- Note that it suffices to find a unitary 
transformation U that transfers I\ x — » I nx and I\ v — > I ny . The same unitary transformation will also 
transfer [I\ Xl Ii y ] — > [Inx, Iny], i-e. I\ z — ► Inz- Therefore, the transfer of the coherent state of spin 1 to 
spin n is equivalent to the transfer of the non-Hermitian operator Zj~ = I\ x —ilx y to I~ = I nx — il n y 
The transfer between such non-Hermitian operators arises naturally in coherent spectroscopy of 
ensembles and constitutes a fundamental step in multidimensional NMR spectroscopy of biological 
macro-molecules . We emphasize again that any unitary transformation that transfers Jj~ to I~ , 
implies that any state of the form po = il + a\Ii x + a-2li y + a^I\ z is transferred to spin n. 

Besides applications in spectroscopy, finding optimal methods to control the dynamics of coupled 
spin networks is of fundamental importance for the practical implementation of quantum information 
processing. In recent years, many innovative proposals have come out to harness the dynamics of 
spins in the liquid |(| and solid state |lj], for the purpose of information processing. Like 
many coherence transfer experiments in multidimensional NMR spectroscopy, these NMR quantum 
computing architectures rely on elaborate sequences of radio frequency (RF) pulses for realizing 
desired effective Hamiltonians. Recent proposals by Yamamoto |fl8[ use a chain of nuclear spins 
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i in the solid state for purpose of computing. A major challenge in such architectures, which is 
also found in various other quantum information devices, is finding efficient ways of making qubits 
interact if they are not directly coupled. A prototype example of this problem is finding efficient 
ways to generate unitary transformations which exchange the states of spins on the two opposite ends 
of a spin chain. Pulse sequences presented in this paper can be used to accomplish such operations 
efficiently. 

Control of Spin Chain Dynamics 

Consider a linear chain of n weakly interacting spin | particles placed in a static external magnetic 
field in the z direction and with Ising type couplings between next neighbors [2C| . In a suit- 
ably chosen (multiple) rotating frame which rotates with each spin at its resonant frequency, the 
Hamiltonian that governs the free evolution of the spin system is given by the coupling Hamiltonian 

n 

Tic = 27T^ Jfc-l,fc I(k-l)zlkz, 
k=2 

where Jk-i >k is the coupling constant between spin k — 1 and k. If the resonance frequencies of the 
spins are well separated, spin k can be selectively excited (addressed) by an appropriate choice of 
the amplitude and phase of the RF field at its resonant frequency. The goal of the pulse designer 
is to make appropriate choice of the control variables comprising of the frequency, amplitude and 
phase of the external RF field to effect a net unitary evolution U(t) which efficiently transfers the 
initial operator A = 1^ to B = I~. For simplicity but without loss of generality (vide infra), we 
assume that all coupling constants in the spin chain are equal, i.e. Jk-i,k = J for 1 < k < n. 

A straight forward way of transferring the operator I-f to I~ is to perform sequential transfers, 
whereby Ij~ is transferred to 7 fc+1 |2^, |j| . Each of these sequential steps takes units of time, 
resulting in a total time of 3 ^j 1 - ) (see Fig.Q). However this is far from optimal. It can be shown 
Q that given any initial operator A k that represents a state involving spin k and spins with label 
less than k, the minimum time required to advance this operator one step up the spin chain is ^j. 
This suggests that an efficient approach to transferring the state to I~ is to prepare an encoded 
state A^~ such that it is possible to go from A^ to A^ , 1 in a time of jj. Furthermore, it is desirable 
that these encoded states be sufficiently localized so that A± can be encoded and decoded in a short 
time. We will refer to such encoded states as effective soliton operators (see end of this section). 

Now we consider the three specific effective soliton operators A^ — c ZI(k-2)xI(k-i)zi ^ky = 
2i(fe-i) x 4z, and A kz = 4J^ k _ 2 ) x I(k-i)yhz, which obey the commutation relations [A a ,A ] = 
ie Q/ 3 7 A 7 , where e Q /3 7 is the Levi-Civita symbol which is 1 (or —1) if {a(3j} is an even (or odd) 
permutation of z} and if two or more of the indices a, /3, 7 are identical. Each individual 
soliton operator A ka is advanced along the spin chain by one unit if the propagator 

U A = cxp{-iAft c } exp{-i|iy (1) 

with F y = Ii y + l2 y + ■■■Iny is applied: 

The propagator U\ can be realized by applying a non-selective 90° pulse (with negligible duration) 
to all spins, followed by the evolution of the spin system under the coupling Hamiltonian TL C for a 
duration A = (2J) _1 . 
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With the help of the soliton operators A fc = Ak x — iAfc a , it is possible to transfer I 1 = I\ x — ili y 
efhciently to /~ = I nx — \I ny : 

17 A- ^ °± A- j- 

1 O ll It 

(n— 3)timcs 

Here, the encoding of 17 as the soliton operator A^ is effected by the propagators 

U\ = exp{-iA74}exp{i-Ii^}exp{-i~I 1:c }, 

and 

U 2 = expj-iAWcJexpj-i-^ia; + h y )}- 

Finally, the decoding of the the soliton operator A~ into I~ is achieved by the propagators Ua and 
U n+ i = exp{i|/„ a; }exp{-iA7i c }exp{i|(/„ 2; - I( n -i) y )}. U\, U 2 , U\, and U n+ i require a period 
A = (2J) _1 each, resulting in the time 

_ n + 1 

7~soliton — 

for the complete transfer from 17 to I~ . 

The flow of soliton operators is summarized in Fig.2B. The panel schematically traces the evo- 
lution of the initial operators I\ Xl I\ y , and I\ z via the soliton operators A^, Ak y , and Ak z in the 
spin chain as a function of time. The operators A/^ , A^j, and A^ 2 represent local correlations of 
spin k — 1 with its neighbors. Under the proposed pulse sequences, these correlations advance one 
step in the spin chain, every T step = jj units of time. Although these operators evolve to other 
operators under the proposed pulse sequences, if the spin system is observed stroboscopically, every 
T st ep units of time, the correlations maintain their shape and are just translated one step up in the 
spin chain. Hence the name effective soliton operators. 



Efficiency of Transfer 



The time T so m on taken by the proposed pulse sequence should be compared with the transfer time 
for conventional pulse sequences which transfer 17 to I~ iffl, 22, 2^|. These pulse sequences require 



n — 1 steps of selective isotropic transfers in which the j th step transfers the operator 17 to IJ +1 . In 
the j step, only spins j and j + 1 are active and the remaining spins in the chain are decoupled. 
This mode of transfer is depicted in panel A of Figure |l]. Each such isotropic transfer step requires 
units of time and therefore the total time is ■ bi the limit of large n, the proposed soliton 

sequences only take | amount of time as compared to state of the art pulse sequences. A comparison 
of the time taken for the coherence transfer by the conventional sequence of selective isotropic pulse 
sequences T conv and the proposed pulse sequences T so uton is shown in the Figure |^ for n < 10. 

The proposed pulse sequences are also compared with the widely used concatenated INEPT 
pulse sequences p4| , which transfer only one component of magnetization along a spin chain, i.e. 
Iix — * Inx- If all couplings are equal to J, the time required for transferring I\ x — * I nx by the 
concatenated INEPT pulse sequences is jj. For large n, this is approximately the same as the 
time required for the new soliton-based pulse sequences. However, the soliton sequences transfer the 
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complete state of spin 1 to spin n, which may result in appreciable gain of signal to noise ratio in 
spectroscopic applications Q. 



The proposed pulse sequences can be used to efficiently exchange the (arbitrary) states of spins at 
the two ends of a spin chain with possible applications to proposed quantum computing architectures 
p8| . This exchange operation between spin 1 and n, in general will not preserve the state of other 
spins on the spin chain and hence donot represent a swap gate between spin 1 and n in the usual 
sense. The proposed soliton sequences can be used to transfer the state of spin 1 to spin n and spin 
n to spin 1 simultaneously in ^j- units of time. This should be compared to the approach where 
states of spin 1 and n are exchanged through a sequence of neighbouring swap operations pl| . Each 
swap operation requires units of time. Therefore the total time required to exchange states of 
spin 1 and n by sequential swapping is at least ^j 1 ^ , which in the limit of large n is three times 
longer than the proposed approach based on soliton operators. 

Finally, we analyze the time taken by the proposed coherence transfer methodology when cou- 
plings are not equal. To produce the effect of the propagator U\ acting on A^, A^, A^ z in equation 
[I], the only terms in the coupling Hamiltonian 7i c that are instrumental are 

2lT{J k _2,k-lI(k-2)zI{k-l)z + Jk-l,kI(k-l) z Ikz + Jk,k+llkzl(k+l)z) ■ 

If J k ,k+i is the smallest of the coupling constants Jk-2.k~i, Jk-i,k> Jk.k+i, i-e. if 

Jk,k+1 = mm (Jfc-2,fc-l, Jfc-l,k) Jk,k+l), 

then it will take Tk — — units of time to produce the propagator 

exp{-i7r(/( fe _ 2)z /( fc _ 1)z +I{k-i)zhz + hzl(k+i)z)}- 

This is achieved by letting the coupling Jk,k+i evolve during Tk while letting Jk-2,k-i and Jk-x,k 
couplings evolve only during Jfc r ~ 2,fc ~ 1 Tfc and J T h ~ 1,k Tk respectively and decoupling these couplings for 
the remaining time. This can be achieved by standard refocusing techniques Q. Therefore the total 
time required for propagation of the soliton A3 to A„ is 

n— 3 ^ 



^-J 2 min(Jfe i fe + i, Jk+i,k+2, Jk+2,k+3) 

Similar arguments yield that the time required for preparation of soliton state from the initial state 

is 1 ^— and finally the time required to reduce the soliton state to the final state is 

2J12 2 mm(j 12 ,j 2 3) 

. 1 + I 

2 min(j„_3,„_i,j„_i, n ) 2 ./„_i, n ' 

The proposed methods for control of spin chain dynamics may have the potential to improve 
the senstivity of multi-dimensional heteronuclear triple resonance experiments, used for example for 
sequential resonance assignments in protein NMR spectroscopy JIo[] . The proposed method of ma- 
nipulating dynamics of spin chains could also reduce decoherence effects in experimental realizations 
of quantum information devices Jig ], Although minimizing the time required to produce desired 
unitary evolutions in a quantum system is expected to reduce dissipation and relaxation effects, 
optimizing pulse sequences by incorporating a realistic relaxation model may further improve the 
sensitivity of experiments in coherent spectroscopy. 
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Figure 1: Panel A shows the flow of coherence and polarization in a spin chain under a sequence 
of selective isotropic mixing periods pj| p2| , p3[ ], (sequential swap operations) each of which can 
be decomposed in three steps of duration A = (2J) _1 with effective Hamiltonians 2irJI( k _ 1 } x Ik x , 

i)zhz, respectively. After the first step the initial operators h 



2t J I (k-i)y hy i and 27rJ7( fe _ 



are transferred to I\ x , 2Ii z l2 X , 21\ y li x respectively. Coherence is transferred in a sequential manner 
where the state of the spin k is transferred to spin k + 1 in M- units of time. The total transfer takes 

units of time. Panel B shows the flow of coherence and polarization under the proposed pulse 
sequence based on effective soliton operators (indicated by grey arrows) . Here, a localized spin wave 
is created which moves one step in the spin chain in every seconds. The total transfer time for 
the proposed pulse sequence is -^j-- For clarity, operators such as I kx are indicated by the letter x 
at position fc. Similarly, bilinear (or trilinear) product operators such as hk x (or Afc 2 ) are indicated 
only by the axis labels x, y, or z at the corresponding spin position, omitting prefactors of 2 (or of 
4) and possible algebraic signs. 
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Figure 2: The figure shows a comparison of r step = t/(ji — 1), the average time required to advance 
by one step in a chain of n coupled spins for pulse sequences which effect full transfer from J{~ to 
Diamonds: conventional sequence of selective isotropic mixing steps between neighboring spins 
|||, H§. Squares: Sequence of optimal indirect SWAP(k, k+2) operations §j| (of duration ^) 
which are followed by a selective isotropic mixing steps between spins (n — 1) and n if n is even. 
Circles: Soliton pulse sequence. 



8 



